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(Dated: January 8, 2013) 

We report the structural and electronic properties of an artificial graphene/Ni(lll) system obtained by the 
intercalation of a monoatomic layer of Ni in graphene/Ir(l 1 1). Upon intercalation, Ni grows epitaxially on 
Ir(lll), resulting in a lattice mismatched graphene/Ni system. By performing Scanning Tunneling Microscopy 
(STM) measurements and Density Functional Theory (DFT) calculations, we show that the intercalated Ni layer 
leads to a pronounced buckling of the graphene film. At the same time an enhanced interaction is measured by 
Angle-Resolved Photo-Emission Spectroscopy (ARPES), showing a clear transition from a nearly-undisturbed 
to a strongly-hybridized graphene 7r-band. A comparison of the intercalation-like graphene system with flat 
graphene on bulk Ni(l 11), and mildly corrugated graphene on Ir(l 1 1), allows to disentangle the two key prop- 
erties which lead to the observed increased interaction, namely lattice matching and electronic interaction. Al- 
though the latter determines the strength of the hybridization, we find an important influence of the local carbon 
configuration resulting from the lattice mismatch. 



I. INTRODUCTION 



Metal supported graphene has received renewed interest as it provides a model-system for studying graphene modifications on 
well-defined large area samples. Recent photoemission and Scanning Tunneling Microscopy (STM) studies have shown that the 
carrier mobility, chirality, and band gap can be tailored by a periodic perturbation potential flj, doping [5], intercalation 0], 
and hybridization with the supporting substrate [5 1. Although there are countless studies on graphene grown on transition metals 
the vastly differing interaction of graphene (G) with transition metals (Me) is not fully understood on a basic level (8l[9]- 
As an important contribution to the interaction stemms from non-local (van-der-Waals like) interactions, the variability of a G 
film with the supporting metal is only partially explained by the so-called d band model [ 10 1 . That Pt and Ir interact more 
weakly with G than Ni or Co is not surprising according to this model, but changes from one metal to the next one in the periodic 
table are expected to be more gradual. Instead, for neighboring elements, like Pd and Pt, and Rh and Ir, the G-Me interaction 
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seems to abruptly switch from strong to weak |5]|6). According to the experimental findings, the G-Me interaction has been 
partitioned into these two main categories, where from an electronic structure point of view, a strong or a weak interaction means 
a perturbed or an almost unperturbed graphene 7r-band at the K-point of the Brillouin zone. On the other hand, the term strong 
appears inappropriate if intended for chemisorption between graphene and the underlaying metal. Indeed, for Ni(l 1 1), which is 
considered one of the metals belonging to the strong category, only a moderate adsorption energy of 67 meV per carbon atom has 
been recently evaluated on the basis of high-level many -body calculations (and up to 160 meV/C using semi -empirical forcefield 
corrections), which still is in the range of typical physisorption systems lfTTHT3l . Therefore, it should be noted that the strong 
interaction mainly implies a strong hybridization between the graphene ir states and the substrate, but is not necessarily reflected 
in the adsorption energies. A related question is a possible correlation between the G-Me lattice mismatch and the strength of 
interaction. Several metals belonging to the strong (Rh, Ru) and weak (Ir, Pt, Cu) categories form moire structures, comprising 
different interactions with the G layer. 

Among several transition metals, Ni(lll) has been most studied as substrate material for G-Me interface, both by theory and 
experiment. The close lattice match between G and Ni allows the growth of a commensurate (lxl) graphene overlayer, with 
carbon atoms at atop and fcc-hollow sites, separated from the substrate by 2.11 A and 2.16 A, respectively lfT4l . Angle-resolved 
photoemission (ARPES) data show a pronounced energy gap at the K-point of the Brillouin zone between ir and ir* 1031 . as a 
result of broken symmetry for the two carbon sublattices accompained by strong hybridization between Ni 3d and graphene ir 
states. 

Here we report a new G-Ni(l 11) system, obtained by the intercalation of a single epitaxial layer of Ni in graphene/Ir( 1 1 1). For 
this system, the lattice mismatch between graphene and the Ni layer is increased. Although the epitaxial growth of Ni on Ir(l 11) 
will lead to additional electronic effects, such as a narrowing of the Ni d band, the local chemical electronic environment is still 
similar enough to allow for a comparision with G/Ni(lll). The intercalation leads to a locally enhanced interaction, resulting 
in a strong corrugation of the graphene layer. We investigated the artificially mismatched G/Ni by STM, Density Functional 
Theory (DFT) including van der Waas contributions (vdW-DF) and ARPES, providing a wide characterization of electronic and 
structural properties. The comparison between G/Ni/Ir(l 1 1) and G/Ni(l 1 1) allows to rank the influence of two important factors, 
lattice mismatch and chemical interaction, affecting the G-Me adsorption mechanism. 

II. METHODS 

The presented studies were performed in two different experimental chambers under identical experimental conditions, allow- 
ing for a reproducible sample preparation. Photoemission experiments were carried out at the VUV-Photoemission beam line 
of the Elettra synchrotron radiation facility (Trieste, Italy), using a Scienta R4000 electron energy analyzer at a base pressure of 
5 x 10~ n mbar. Angle-resolved photoemission spectra were collected at room temperature (RT) using a photon energy of 80 
eV, with total energy resolution of 100 meV and angular resolution of 0.1°. 

STM experiments were carried in ultra-high vacuum (UHV) system (base pressure 5 x 10 -11 mbar) equipped with an Omicron 
variable temperature scanning tunneling microscope. All STM measurements were performed in the constant-current-mode at 
RT using electrochemically etched polycrystalline tungsten tips cleaned in UHV by flash-annealing. The sign of the bias voltage 
corresponds to the voltage applied to the sample. Tunneling current and voltage are labeled It and Ut, respectively. 
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FIG. 1 : (color online) (a) Topographic STM overview showing the morphology of graphene with a partially intercalated Ni sub-monolayer. 
Ni accumulates at step edges (B area) showing increased moire corrugation in STM as compared to G/Ir(l 1 1) (A and C areas) (70 x 70 nm 2 ; 
Ut = 0.65 V; It = 1.21 nA). Corresponding LEED image in the inset, (b) Areas with Ni intercalated underneath graphene (B areas) show 
reduced mean apparent height in the line profiles and the histogram. The histogram shows the frequency of apparent height values appearing 
in the magnification depicted in (c) (black curve) and within areas on terrace A, B or C (yellow, orange and brown curves, respectively), (c) 
Magnification of the dotted square in (a) (46 x 8.6 nm 2 ). 

Experimentally, the G/Ir(lll) system was prepared by the procedure described in JT|. Intercalation of Ni underneath a 
graphene layer was performed via annealing of the pre-deposited film in the temperature range of 670-800 K. Starting from the 
sub-monolayer regime, the Ni coverage was estimated on bare Ir(l 1 1) by measuring the intensity ratio of Ni-3p and Ir-4f core 
levels. 

The spin-polarized DFT calculations were performed with the Vienna Ab-initio Simulation Package (VASP) fT6l[T7l . using 
PAW potentials lfl8l \l9l and an energy-cutoff of 400 e V. As GGA exchange-correlation functionals tend to severely underbind 
the adsorption of graphene on Ni(lll) [11], the calculations were performed using van der Waals DFT (vdW-DF) with the 
opt86b functional ||20l [2TI . In the calculations, a (10 x 10) graphene sheet was adsorbed on a (9 x 9) Ir(lll) substrate with 
a lattice constant of 2.735 A, consisting of a three layer slab and an additional intercalated epitaxial Ni layer. A T-centered 
3x3x1 k-point mesh was used to relax the structures keeping only the two bottom-most layers fixed. The C Is core level 
shifts were calculated in the initial state approximation. For the graphical visualisation, the resulting total core level spectra are 
displayed as a sum over Gaussian functions with a standard deviation of 0.25. The STM simulations were performed using the 
Tersoff-Hamann approximation ll22l using the integrated charge density between Ep and Ep + 0.2 eV. 



In STM, the G/Ir(lll) surface displays large fully graphene covered terraces with several hundreds of nanometers width and 
straight steps following the direction of the graphene moire, consisting of distinct fcc-hcp carbon configurations and virtually 
indistinguishable top-hollow (top-fcc and top-hep) sites. G/Ir(lll) was imaged here in the dark-atop-contrast [23], where ele- 
vated fcc-hcp regions appear as black depressions in middle of bright rings in STM topographies. Upon Ni intercalation straight 
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FIG. 2: (color online) (a) Structural model for a single layer of graphene on Ni/Ir(l 11). The color coding indicates the height of the 
corrugation Ah in the graphene layer, (b) Comparison of the corrugation in the optimized structure of graphene/Ni/Ir(l 1 1) (upper panel) and 
graphene/Ir(l 1 1) (lower panel), (c) Simulated STM image for the states between Ep and 0.2 eV. (d) Atomically resolved STM topography of 
graphene/Ni/Ir(lll) (8x8 nm 2 ; U T = 50.0 mV; I T = 35.0 nA). Theoretical (e) and experimental (f) Cls core level shifts of G/Ni/Ir(lll) 
compared to G/Ir(l 1 1). 



terraces become disrupted and irregularly extended by areas with inverted moire contrast (white protrusions on dark background) 
as compared to pristine G/Ir(l 1 1) steps. Figure la depicts a detailed STM topography with graphene covering Ir (A, C) including 
an intermediate area with inverted contrast (B) and two remaining Ni clusters (in white) on top of graphene with height in the 
nanometer range. From the morphology of the sample after intercalation it becomes clear that areas A and C display pristine 
graphene on adjacent Ir substrate levels, whereas area B corresponds to graphene on a Ni-intercalated region. To shed more light 
on the overall as well as site specific graphene-substrate interaction we analyzed the area depicted in Figure lc and evaluated 
line profiles across the terraces and histograms showing the distribution of apparent height values (Figure lb). 

Line profile 1 crosses the fcc-hcp sites of G/Ir(lll) on terrace A, which appear as dark depressions, and continues across 
terrace B where bright protrusions now occupy the former fcc-hcp sites. Height profile 2 crosses the bright protrusions of terrace 
B. A large peak-to-valley height variation of 0.6 A is measured on terrace B as compared to the much shallower peak-to-valley 
corrugation of 0.25 A of G/Ir(l 1 1) on terraces A and C. In the histogram in Figure lb the frequency of apparent height values is 
shown for equally sized areas on terrace A, B and C, respectively, as well as for the complete area in Figure lc. For G/Ir(l 1 1) on 
terrace A and C the distribution is narrow (0.45 A peak width) and features a distribution maximum reflecting the top-hollow sites 
and a distinct shoulder at lower apparent height corresponding to the fcc-hcp regions. For Ni intercalated graphene on terrace 
B the distribution is much wider (0.8 A peak width) with a maximum 0.6 A below the maximum of G/Ir(lll) and a shoulder 
extending far into the G/Ir(l 11) region. Assuming the intercalated Ni atoms arrange pseudomorphically on the Ir( 111) surface 
with comparable interplane distance, the measured distance between equivalent points on the terraces A and B reflects to a large 
extent the difference in the graphene-metal distance (see also the discussion of the DFT results), and leads to the intriguing result 
that graphene on intercalated Ni is in average 0.6 A closer to the topmost substrate layer compared to the pristine G/Ir(lll). 
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The almost unaffected continuation of the graphene moire on intercalated Ni - albeit with increased corrugation - justifies the 
assumption of pseudomorphic arrangement of the intercalated Ni atoms. Moreover, the experimental data do not show any 
indication for the formation of a surface alloy. 

The DFT calculations allow to investigate the structural changes in the (10 x 10) graphene sheet adsorbed on a (9 x 9) 
Ni/Ir(l 1 1) substrate [1 ML Ni pseudomorphically arranged on Ir(l 1 1)]. It should be noted that the strength of the hybridization 
is closely related to the minimal graphene-substrate distances ifTTIl . Consequently the structural analysis allows to deconvolute 
the effects of the lattice mismatch between the graphene sheet and the substrate, and the chemical properties of the interface. 
Figures 2(a-b) display the geometry of the graphene sheet after the relaxation. The model shows that the intercalation of the Ni 
layer leads to a pronounced corrugation of Ah = 1.51 A in the graphene layer, significantly larger than for G/Ir(lll) (Figure 
2b). Yet it should be noted that more than 70% of the carbon atoms in the graphene layer are adsorbed at a close distance of 
about 2.0 - 2.2 A from the Ni layer. The comparison with the inter-plane distance of G on Ni(lll) (2.1 A) therefore hints at a 
similar binding of graphene to the substrate, despite the large experimental strain of roughly 9% due to pseudomorphic growth 
of the Ni lattice in G/l ML Ni/Ir(l 1 1). Therefore, the G-Ni bonding seems to be mainly affected by the electronic contribution, 
which also drives the strong corrugation of the graphene layer in the mismatched structure. 

In the flat regions, the magnetic moment of the Ni atoms is completely quenched by the interaction with the graphene sheet, 
while the Ni atoms under the graphene bubbles yield a small magnetic moment (<0.4 /is). We find that no magnetic moment 
is induced in the graphene sheet that can be due to the small net magnetic moment of the underlying Ni film, opposite to the 
graphene/Ni(l 1 1) system showing an induced magnetic moment of carbon atoms 11241 . 

Seen from an atomistic point of view, the close adsorption configuration of the graphene layer is reached not only for the 
top-fcc configuration preferred on Ni(l 1 1), but also for the adsorption in a bridge-like configuration. Both configurations yield a 
close adsorption distance reaching values as low as 1 .94 A. In contrast, the weak interaction in the fcc-hcp sites (green regions in 
Figure 2a) leads to the formation of local protrusions, with a maximal distance of 3.45 A to the Ni layer common for physisorbed 
graphene. Nevertheless, this distance is still smaller than the calculated maximal (vdW-DF) distance 3.7 A (Ah = 0.37 A) for 
the adsorption of G on the bare Ir( 111) surface (Figure 2b). 

A direct comparison of the obtained STM data (Figure 2d) and a simulated STM image (Figure 2c) reflects the structure of the 
adsorbed graphene sheet: the elevated fcc-hcp regions appear brightest, while the low-lying areas with top-hollow configuration 
appear as a dark background. In agreement with the structure recently reported for G/Ru(0001)[25 1 and G/Rh(l 1 1) [26-28), the 
regions where the graphene sheet is adsorbed in a local bridge configuration is the area of the smallest distance to the surface. 
These areas appear as faintly visible depressions in STM topographies. A peak-to-valley corrugation of up to 1 A fits well the 
corrugation of 1.3 A in the simulated image. On the atomic level, G/Ni/Ir(l 1 1) shows rings of carbon atoms everywhere within 
the moire supercell in Figure Id, however with the strongly bound areas a difference in the intensity between the neighboring 
atoms is observed indicating a broken sublattice symmetry. 

Previous studies demonstrate that the corrugation and hybridization of the graphene layer with the metallic substrate is strongly 
reflected in the Cls line-shape ll29l[30l . Figure 2f shows the Cls core level (CL) taken at 445 eV of G/Ir(lll) and its evolution 
during the intercalation of about one third of monolayer and one monolayer of Ni atoms. According to the existing literature 
H9][31), in the G/Ir(l 11) system the Cls binding energy is found at (284.10±0.20)eV. After the intercalation of 0.33 ML of Ni 
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FIG. 3: (color online) (a-d) ARPES dispersions along TK as a function of the amount of Ni atoms intercalated underneath G/Ir(l 1 1). (a) and 
(d) show extreme cases of ML and thick Ni, respectively, (e) Carbon projected band structure of a (lxl) model system of graphene/Ni/Ir(l 1 1). 
The band structure was evaluated at a G-Ni distance of 2.0 A (brown-dark gray dots) and 3.4 A (green-light gray dots). 

atoms a second component at higher binding energy is seen. The relative intensities of the two components is fully inverted for 
1ML of Ni atoms intercalated. The main peak centered at 284.90 (±0.20) eV is close to the value found for graphene grown 
on bare Ni(lll) l32l , where a single peak at 284.7 (±0.18) eV was measured, with an intrinsic line width of 216 meV. In our 
system, the Cls line-shape exhibits a total width of about 840 meV and a strong asymmetry towards lower binding energy, likely 
convoluting different components. This is confirmed by the calculated core level states of G/Ni/Ir(l 1 1): although the calculations 
predict only a minor CL shift for the carbon atoms in the elevated regions of the moire pattern, the strongly interacting C atoms 
in the lower regions are dominant in the convolution considering all contributions (Figure 2e) and thus do not exhibit a double 
Cls peak as observed for G/Re(0001), G/Rh(l 1 1) or G/Ru(0001) I29ll30l . 

To shed more light on the overall interaction of graphene with the mismatched Ni layers, we have performed ARPES mea- 
surements mainly looking at the graphene 7r-band. Figures 3(a-d) show ARPES maps of the electronic band dispersion of (a) 
G/Ir(lll); (b) G/0.33MLNi/Ir(lll); (c) G/lMLNi/Ir(lll); (d) G/fhickNi/Ir(lll), all taken along the TK direction. The last 
measurement was collected after the intercalation of several monolayers of Ni. The tt state of G/Ir(lll) (Figure 3a) exhibits a 
minimum at the T-point at 8.30 eV, and approaches the Fermi level at the K-point at 70meV, where the planar a state reaches 
11.35 eV of binding energy (TJ. According to the literature [1,3 29] [33] [34), replicas bands of both tt and a states due to the 
moire superpotential are seen close to the K-point. After the intercalation of about 1/3 ML of Ni atoms (Figure |3j)), new tt 
and a states appear at higher binding energy together with the d states of Ni, weakly dispersing in the range 0-2 eV below the 
Fermi level. The co-existence of double tt and a states reflects a non-homogenous surface with clean areas of G/Ir(lll) and 
patches where Ni atoms are in between. When a full monolayer of Ni atoms is intercalated via annealing (Figure 3c), the new 
band structure evolves clearly. The 7r state now exhibits a minimum at the r -point at 10.03 eV and reaches a maximum at about 
2.16 eV at the K-point, where it merges with the d states of Ni. The a state is also shifted to a higher binding energy compared 
to G/Ir(lll), with a minimum at the K-point at about 12.48 eV. When several monolayers of Ni (above 5) are intercalated via 
annealing (Figure 3d), the electronic states of Ir are no longer detected, and the band structure reflects that of G/Ni(lll) lfI31 . 
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Notably, while the a state is unaffected by the number of Ni layers at the K-point, the tt state exhibits a clear shift towards lower 
binding energies (up to 2. 16 eV for 1 ML intercalated Ni) with respect to G grown on multilayers (Figure 3d) or bulk nickel [15], 
where the tt band maximum is found at about 2.65 eV. This finding is related to the differences in width of the Ni-d states of a 
single intercalated Ni layer compared to a surface of bulk nickel: narrowing of the Ni 3d band upon the decreasing of the Ni layer 
thickness. Taking into account that the position of the graphene-derived tt band at the T point is the same for both thick and thin 
(1 ML) intercalated Ni layers, we can conclude that the energy shift of tt band with respect to free-standing graphene is purely 
defined by the charge transfer between Ni and C atoms at the closest distance through the donation/back-donation mechanism 
(9). At the same time the presence as well as the width of the band gap between tt and tt* graphene-derived states is determined 
by the broken symmetry for two carbon atoms in the graphene unit cell in this system accompanied by a strong hybridization 
between Ni 3d and graphene tt states. 

For a comparison with experiment, we also evaluated the theoretical band structure. In order to avoid the back-folding 
induced by the larger supercell, we have mimicked the local interactions by calculating the electronic structure for a smaller 
(lxl) model of G/Ni/Ir(lll) in a top-fcc configuration at an average distance of the flat regions (2.0 A) and at the maximal 
height of the bubbles (3.4 A). The resulting band structure is shown in Figure 3e. Although the effect of the lattice mismatch 
between the graphene sheet and the substrate is lost in this smaller model, the calculations clearly show that the interactions 
at the elevated regions of the bubbles (green-light gray dots in Figure 3e) are rather weak, resulting in a nearly unperturbed 
graphene band structure. On the other hand, a much stronger interaction can be expected for the dominant flat regions in the 
vicinity of the surface, leading to a large splitting of the tt band at the Dirac point (brown-dark gray dots in Figure 3e). These 
findings agree with the experimentally observed opening of a band gap in the ARPES data. Furthermore, no tt band splitting 
is observed experimentally in our system, due to the metallic nature of graphene 11351 . in contrast to the electronic behavior of 
h-BN grown on selected transition metals l36l . where the dielectric nature of the overlayer allows to observe double a and tt 
states corresponding to upper and lower regions. 



IV. CONCLUSION 



In conclusion, we have shown that the adsorption of graphene on epitaxial layers allows to study the influence of the lattice 
mismatch between the graphene layer and the support, while keeping the chemical environment similar. For graphene on 
Ni/Ir(lll), we find that the interaction is locally strongly enhanced for specific adsorption configurations. Consequently, in 
contrast to G/Ir(lll), the moire structure of G/Ni/Ir(lll) exhibits a strong corrugation, with a modulation of about 1.5 A and 
a minimum G-Me distance slightly smaller than 2 A . The graphene band structure probed by ARPES shows a clear transition 
from a nearly-free standing to a strongly-hybridized character of the tt band, in analogy with graphene grown on bulk nickel. 
The hybridization between Ni d states and graphene tt states is directly related to the strongly interacting top-hollow and bridge 
configurations in the lower parts of the moire mesh, at a distance of about 2.0 - 2.2 A from the Ni layer. In contrast, the interaction 
is significantly weaker for other regions (fcc-hcp configurations) of the moire mesh, where only a van der Waals like binding is 
observed. Therefore we can identify the role of two important contributions to the adsorption: while the electronic interaction 
dominates in the strongly interacting regions, the lattice mismatch between graphene and the metal support is decisive for the 
ratio between strongly and weakly interacting regions. We expect that this interplay is one of the key features for mismatched 
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graphene-metal interactions. 
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